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Abstract — This paper presents Matlab script based 
simulation model for a PV cell in order to allow estimate the 
electrical behavior of the cell with respect changes on 
environmental conditions of temperature and irradiance. An 
accurate PV module electrical model is presented based on the 
Shockley diode equation. The general model was implemented 
using Matlab script, and effect of irradiance and temperature as 
variable parameters and effect on outputs the I-V characteristic 
is studied. A particular typical 180W solar panel was used for 
model evaluation, and results was compare with points taken 
directly from the manufacturer’s published curves and show 
good correspondence to the model. 

Index Terms — Photovoltaic, photocurrent, photovoltaic 
module, insolation. 


I. Introduction 

Nonconventional energy resources will be an increasingly 
important part of power generation in present case. 
Regardless of the intermittency of sunlight, solar energy is 
widely available and completely free of cost. Recently, 
photovoltaic array system is likely recognized and widely 
utilized as important source in electric power applications [1]. 
Besides assisting in the reduction of the emission of 
greenhouse gases, they add the much- needed flexibility to the 
energy resource mix by decreasing the dependence on fossil 
fuels. Photovoltaic (PV) systems produce DC electricity when 
sunlight shines on the PV array, without any emissions. The 
DC power is converted to AC power with an inverter and can 
be used to power local loads or fed back to the utility. 

The PV application can be grouped depending the scheme of 
interaction with utility grid, grid connected, stand alone, and 
hybrid. PV systems consist of a PV generator (cell, module, 
array), energy storage devices (such as batteries), AC and DC 
consumers and elements for power conditioning. This paper is 
about a model for modeling and simulation of PV module 
based on Shockley diode equation. 

II. IPV GENERATOR 

Solar cell is basically a p-n junction fabricated in a thin 
wafer or layer of semiconductor. The electromagnetic 
radiation of solar energy can be directly converted electricity 
through photovoltaic effect. Being exposed to the sunlight, 
photons with energy greater than the band-gap energy of the 
semiconductor are absorbed and create some electron-hole 
pairs proportional to the incident irradiation. Under the 
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influence of the internal electric fields of the p-n junction, 
these carriers are swept apart and create a photocurrent which 
is directly proportional to solar insolation. PV system 
naturally exhibits a nonlinear I-V and P-V characteristics 
which vary with the radiant intensity and cell temperature. 


III. MODELLING OF SOLAR CELL 

Simple circuit having current source and parallel diode is 
used for modeling. The shunt resistance R sh is neglected. A 
single shunt diode was used with the diode quality factor set to 
achieve the best curve match[2].This model is a simplified 
version of the two diode model. The circuit diagram for the 
solar cell is shown in Figure 1. 
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• Shunt resistance R sh , in parallel with the diode, this 
corresponds to the leakage current to the ground and it is 
neglected 
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I =Diode saturation current, 
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I L =Photocurrent 

Rs=Series resistance 
n= diode quality factor 

Considering importance effect of temp and solar radiation can 
be introduced to the model by adding in turn 

• Temperature dependence of the diode saturation current I 0 . 

• Temperature dependence of the photo current I L . 

• Series resistance R s , which gives a more accurate shape 
between the maximum power point and the open circuit 
voltage. This represents the internal losses due to the current 
flow 

The model included temperature dependence of the 
photocurrent I L and the saturation current of the diode I 0 [ 3]. 
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A series resistance Rs was included; which represents the 
resistance inside each cell in the connection between cells. 
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Fig.3 P-V Characteristics of PV module at 25°C 

Fill factor is the ratio of the maximum power that can be 

delivered to the load and the product of Isc and Voc: 
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The fill factor is a measure of the real I-V characteristic. Its 
value is higher than 0.7 for good cells. The fill factor 
diminishes as the cell temperature is increased. The open 
circuit voltage increases logarithmically with the ambient 
irradiation, while the short circuit current is a linear function 
of the ambient irradiation. The dominant effect with 
increasing cell’s temperature is the linear decrease of the open 
circuit voltage, the cell being thus less efficient. The short 
circuit current slightly increases with the cell temperature. 


IV. I-V AND P-V CHARACTERISTICS OF SOLAR 

MODULE 

A typical I-V characteristic of the solar cell for certain 
ambient irradiation G and a certain fixed cell temperature T aC 
is obtained using script fig.2. For a resistive load, the load 
characteristic is a straight line with scope IPs/-MR. It should 
be pointed out that the power delivered to the load depends on 
the value of the resistance only. 



Fig.2 A typical current-voltage characteristics of PV module 


Maximum power point is the operating point A ( P max 7 max ) 
in Fig 3, at which the power dissipated in the resistive load is 

p =VI 


maximum: 


max 


max max 


V. EFFECT OF TEMPERATURE AND IRRADIATION 
ON OUTPUT OF PV MODULE 

The influence of the ambient irradiation G and the cell 
temperature T aC on the cell characteristics , can be obtained 
from the model equations. The PV cell photocurrent 7 L (A) is 
directly proportional to solar irradiance G(W/m2). When the 
solar cell is short circuited, negligible current flows in diode. 
Hence the proportionally constant in (3) is set to the rated 
short circuit current I S c at is delivered under rated irradiation. 
Solar intensities are commonly normalized with respect to full 
solar radiation at sea level with average humidity and aerosol 
particulate concentration (1 Sun = 1000Watt/m2). Though 
somewhat contrary to intuition, PV cell performance does not 
degrade significantly between full sun and cloudy conditions. 
The power output decreases lineally with incident solar 
energy, but efficiency is nearly flat over the region of concern. 
The relationship between the photo-current and temperature 
is linear (eqn. 2) and is deduced by noting the change of 
photo-current with the change of temperature (eqn.4). When 
the cell is not illuminated, the relationship between the cell’s 
terminal voltage and current is given by the Shockley 
equation. When the cell is open circuited and illuminated, the 
photo-current flows entirely in the diode. The I-V curve is 
offset from the origin by the photo generated current I L (eqn. 
1). The value of the saturation current I 0 at 25°C is calculated 
using the open circuit voltage and short circuit current at this 
temperature (eqn. 6). An estimate must be made of the 
unknown “ideality factor”. The value of n is between 1 and 2, 
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being near one at high currents, rising towards two at low 
currents. A value of 1.3 is consider as typical in normal 
operation, and may be used initially, until a more accurate 
value is estimated later through curve fitting. The relationship 
of I 0 to temperature is complex, but fortunately contains no 
variables requiring evaluation (eqn.5). The series resistance 
of the panel has a large impact on the slope of the I-V curve at 

V= V oc .[4] 

The TP 180 PV module was chosen for modeling, due is 
well-suited to traditional applications of photovoltaic’s. The 
TP 180 module provides 180 watt of nominal maximum 
power, and has 72 series connected polycrystalline silicon 
cells. The key specifications are shown in Table 1. The model 
of the PV module was implemented using a Matlab program. 
The model parameters are evaluated during execution using 
the equations listed on the previous section. The program, 
calculate the current /, using typical electrical parameter of 
the module ( I sc , Voc ), and the variables Voltage, Irradiation 
(G), and Temperature ( TaC ). The TP180, a typical 180W PV 
module, was chosen for modeling. The module has 72 series 
connected polycrystalline cells. The key specifications are 
shown in table 1 [5]. 

TABLE 1 

TYPICAL ELECTRICAL PARAMETRES OF TP 180 AT STP 


Power Output Pmax(W) 

180 

Power tolerance Pmax(W) 

[-0&+5] 

Module efficiency r)(%) 

14.4 

Voltage at P M ax (V) 

35.8 

Current at Pmax(A) 

5.03 

Open circuit voltage V 0 c (V) 

43.6 

Short circuit current I S c (A) 

5.48 


The model was evaluated using Matlab. The model 
parameters are evaluated during execution using the 
equations listed in the previous section using the above data 
points contained in the script. The current I a is then evaluated 
using these parameters, and the variables Voltage, Irradiation, 
and Temperature. If one of the input variables is a vector, the 
output variable (current) is also a vector. The inclusion of a 
series resistance in the model makes the solution for current a 
recurrent equation (refer to eqn. 1). A simple iterative 
technique initially tried only converged for positive currents. 
The Newton Raphson method used converges much more 
rapidly, and for both positive and negative currents. A listing 
of the Matlab script which implements the equations shown is 
given in Figure 4. 

function I = TP_180(Va,G,TaC) 

% Define constants 

k= 1.381e-23; % Boltzmann’s constant 
q = 1.602e-19; % Electron charge 

% curve fitting of I-V character (Use data for 1000W/m A 2) 
n = 1.62; % Diode ideality factor (n), 

Eg = 1.12; % Band gap energy; 

Ns = 72; % Series connected cells (TP180, 72 cells) 


TrK = 298; % Reference temperature (25C) in Kelvin 
Voc_TrK = 43.6 /Ns; % Voc (open circuit voltage per cell) @ 
temp TrK 

Isc_TrK = 5.48; % Isc (short circuit current per cell) @ temp 
TrK 

a = 0.65e-3; % Temperature coefficient of Isc (0.065%/C) 

% Define variables 

TaK = 273 + TaC; % Module temperature in Kelvin 
Vc = Va / Ns; % Cell voltage 
% Calculate short-circuit current for TaK 
Isc = Isc_TrK * (1 + (a * (TaK - TrK))); 

% Calculate photon generated current @ given irradiance 
I L = G * Isc; 

% Define thermal potential (Vt) at temp TrK 
Vt_TrK = n * k * TrK / q; 

% Define b = Eg * q/(n*k); 
b = Eg * q /(n * k); 

% Calculate reverse saturation current for given temperature 
Ir_TrK = Isc_TrK / (exp(Voc_TrK / Vt_TrK) -1); 

Ir = Ir_TrK * (TaK / TrK) A (3/n) * exp(-b * (1 / TaK -1 / 
TrK)); 

% Calculate series resistance per cell (Rs = 5.1mOhm) 
dVdI_Voc = -1.0/Ns; % Take dV/dl @ Voc from I-V curve of 
datasheet 

Xv = Ir_TrK / Vt_TrK * exp(Voc_TrK / Vt_TrK); 

Rs = - dVdI_Voc - 1/Xv; 

% Define thermal potential (Vt) at temp Ta 
Vt_Ta = n * k * TaK / q; 

% I = I L - Ir * (exp((Vc + I * Rs) / Vt_Ta) -1) 

% f(l) = I L -1- Ir * ( exp((Vc + I * Rs) / Vt_Ta) -1) = 0 
% Solve for Ia by Newton's method: I 2 = Ii - f(Ii)/f (Ii) 
Ia=zeros(size(Vc)); % Initialize I with zeros 
% Perform 5 iterations 
for j=l:5; 

I = I - (I L -1 - Ir .* ( exp((Vc + I .* Rs)./ Vt_Ta) -1))... 

./ (-1 - Ir * (Rs ./ Vt_Ta) .* exp((Vc + I .* Rs)./ Vt_Ta)); 
end 

Fig. 4 Matlab scrip file for the PV module model 

VI. RESULTS OF MATLAB™ PV MODULE MODEL 

The output of the Matlab function is shown first for various 
temperatures (Fig.5) and then for various irradiation levels 
(Fig. 6). These are points taken directly from the 
manufacturer’s published curves, and show good 
correspondence to the model[6]. 



Fig.5 Effect of temperature in output of PV module 
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Fig.6 Effect of solar isolation on PV module 


VII. CONCLUSIONS 

In the present work a Matlab script file for simulation of the 
operation of a photovoltaic system. Software includes the 
models of the solar panel. A generalized PV model which is 
representative of the all PV cell, module, and array has been 
developed with Matlab script file and been verified with a PV 
cell and a commercial module. The proposed model takes 
sunlight irradiance and cell temperature as input parameters 
and outputs the I-V and P-V characteristics under various 
conditions. Such a generalized PV model is easy to be used 
for the implementation on Matlab/Simulink modeling and 
simulation platform. Especially, in the context of the Simulink 
tool, there is now a generalized PV model which can be used 
for the model and analysis in the field of solar PV power 
conversion system. 
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